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Two component Bose gas in a double well potential with repulsive interactions may undergo a 
phase separation transition if the inter-species interactions outweigh the intra-species ones. We 
analyze the transition in the strong interaction limit within the two-mode approximation. Numbers 
of particles in each potential well are equal and constant. However, at the transition point, the 
ground state of the system reveals huge fluctuations of numbers of particles belonging to the different 
gas components. That is, probability for observation of any mixture of particles in each potential 
well becomes uniform. 
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I. INTRODUCTION 

Ultra-cold dilute gases of bosonic atoms constitute per- 
fect systems for experimental and theoretical investiga- 
tions of various phenomena of quantum many body prob- 
lems [l| . From the viewpoint of quantum computing and 
interfcrometry an especially relevant subject is quantum 
fluctuations 0, II] • 

Most experimental studies of fluctuations concentrated 
on systems of cold atoms in double well Q and opti- 
cal lattice potentials [sj. In the former system squeezed 
states were predicted and produced, with particle number 
fluctuations (i.e. uncertainties of populations of the po- 
tential wells) turning from poissonian to sub-poissonian 
. The latter system reveals a superfluid to Mott insu- 
lator transition 



with enhanced phase fluctuations 
but with decreasing particle number fluctuations. 

In the present paper we focus on a system where the 
total particle number is fixed but occupation of certain 
single particle states reveals considerable quantum fiuc- 
tuations. We are interested in a system where the mean 
field theory predicts symmetry breaking [l^ - [T^ and the 
symmetry broken solutions are degenerated and form a 
Hilbert subspace parameterized by a continuous param- 
eter. If the occupation of single particle states varies 
a lot as we move in the degenerate subspace than huge 
particle number fluctuations can be expected in the ex- 
act quantum many body eigenstates. Degenerate sub- 
space parameterized by a continuous parameter appears 
in spin-1 Bose gas with an anti- ferromagnetic interac- 
tion [l5| or in scalar condensates with solitonic solutions 
p^5l - [l7| . Attractive single component Bose gas in a sym- 
metric double well potential reveals also huge particle 
number fluctuations but it constitutes a slightly different 
example [H, [l^. There the degeneracy is smaU, i.e. the 
degenerate subspace is two dimensional, and the particle 
number fluctuations correspond to random localization 
of all particles in one of the potential wells in different 
experimental realizations. In all these examples the cor- 



rect mean fleld theory reduces to the Gross-Pitaevskii 
equations [l|. In the present paper we consider a Bose 
gas system where the Gross-Pitaevskii equation is not a 
correct mean-fleld description, that is, a two-component 
Bose gas in a double- well potential in the strong interac- 
tion limit. 

In Sec. |TI] we present a theoretical model for a two- 
component Bose gas in a double-well potential. In 
Sec. nil Al we derive the effective Hamiltonian using sec- 
ond order perturbation theory valid in the strong inter- 
action limit. In Sec. HUB] we analyze its mean fleld (clas- 
sical) limit and identify phase transition region. It turns 
out that the mean-field solutions reveal continuous de- 
generacy at the transition point. We deduce the exact 
ground state of the system in Sec. IIII CI and show that 
the particle number fluctuations are indeed huge at the 
critical point. In Sec. IIII Dl we estimate the range of 
parameters where the predicted fluctuations can be ob- 
served and in Sec. IIVI the results presented in the paper 
are summarized. 



II. THE MODEL 

The Hamiltonian of a two component Bose gas in a 
symmetric double well potential, in the tight binding ap- 
proximation, takes the form of the Bose-Hubbard model 
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where we have assumed that intra-species interactions 
are the same in both gas components and they are char- 
acterized by a coupling constant U . The parameter Uab 
is a coupling constant that describes inter-species inter- 
actions while J stands for the tunneling rate between the 
two potential wells. We assume also that numbers of par- 
ticles of each component are equal to 2N . Such a choice 
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of the system parameters allows us to perform fully ana- 
lytical calculations. Analysis of a general case is beyond 
the scope of the present paper. The Hamiltonian ([T]) can 
be transformed to 
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where Us = U + Uab, Ud = U — Uab and constant terms 
have been omitted. In the following we consider Us as 
the unit of energy. 
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III. PERTURBATION APPROACH 

A. The second order effective Hamiltonian 

We are interested in the strong interaction limit. 
Therefore, the tunneling part of the Hamiltonian will be 
considered as a small perturbation. For J = the system 
Hamiltonian has exact cigenstatcs 

\N + na,N -na)\N + nb,N -m), (3) 

where N + Ua, N ~ Ha refer to numbers of particles of 
the component a in the first and the second potential 
well, respectively, and similarly for the component b. The 
energies of such states (remember that Us is the unit of 
energy) arc 

E=^{na+nbf + ^{na'nbf. (4) 

Switching to variables Hs — Ua + Ub, rid — Ua ~ nb we 
obtain eigencnergies in a very simple form 

E-lnl + ^nl (5) 

If we assume that the parameters satisfy the condition 

1 » \Ud\N^, (6) 

then manifolds with different values of \ns \ are separated 
on the energy scale (see Fig. [T|). The lowest energy man- 
ifold is related to ris = and states within each manifold 
are labelled by different values of Ud- 

Matrix elements of the tunneling part of the Hamilto- 
nian arc zero between states of the same manifold. How- 
ever, this part of the Hamiltonian introduces couplings 
between different manifolds. In an effective Hamiltonian 
that describes the lowest manifold of the system the ef- 
fect of the coupling can be included via the second order 
perturbation theory. A compact form of the effective 
Hamiltonian may be obtained if we introduce spin oper- 
ators 



FIG. 1: (Color online) Energy levels ([5]) versus Ud for a num- 
ber of particles of each component equal 2N — 100. Black 
lines; the lowest energy manifold, i.e. Us = 0, red lines: the 
manifold corresponding to \ns \ = 1. 

Sjy = -^{aybj~b]aj), 

= i(a]a,-6]5,). (7) 

States belonging to the lowest manifold {ug = 0) can be 
written in the Fock basis ^ as 

N 

1^) = ^{n)\N + n,N -n)\N -71, N + n). (8) 

n=-N 

The Fock states \N + n,N - n)\N - n, N + n) arc the 

cigenstatcs of the S^-, Siz and operators with the cor- 
responding eigenvalues N{N n and — n, respectively, 
so the Us — Q manifold can be specified by: 

5i. + ^2. -0, S2=iV(iV + l). (9) 

In the second order in J the effective Hamiltonian that 
describes the lowest manifold reads (20l - [2^ 

H2 = -2J^^i-^2 + Ud{Sl,_+Sl,). (10) 

The above Hamiltonian together with the condition © 
defines our problem where each potential well is associ- 
ated with an angular momentum operator and there is 
interaction between such subsystems due to tunneling of 
atoms. Note that eigenstates of the Hamiltonian (fTO|) 
depend on two parameters only, i.e. Ud/{2J'^) and N. 



B. Classical limit 

Let us analyze the Hamiltonian ([TO)) [in the manifold 
defined in in the classical limit by substituting the 
spin operators by classical angular momentum compo- 
nents. The condition ^ implies that Siz = —S2z- We 
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are interested in the ground state of the system, 
value of the tunneling part of the Hamiltonian 



The 



— 2,/^Si • S2 — —2j'^{SlxS2x + SlyS2y + 5*12522), (11) 

is minimal for Si^ = S2X, Siy = For Ud/ (2 J^) + 1 < 
the ground state corresponds to \Sjz\ = N which can 
be related to Ud = ±2N (i.e. phase separation occurs 
where different gas components occupy different poten- 
tial wells). For Ud/{2J'^) + 1 > the z-components 
of the angular momenta disappear in the ground state 
{Sjz = 0) which corresponds to = (i.e. equal mix- 
ture of both components in each potential well). When 
Ud/{2J'^) -f 1 wc deal with the transition point where 
1 5^2! can be arbitrary provided 5*12 + 5*22 = 0. Then 
all values of Ud are equally probable (any mixture of 
both components in each potential well is equally likely). 
The transition between the phase separation and miscible 
regimes is discontinuous. 



C. Quantum ground state 

The analysis of the classical limit suggests that huge 
particle number fluctuations can be expected in the quan- 
tum ground state of the system at the transition point. 
Let us switch now to quantum analysis of the Hamilto- 
nian ([T0|) . For Ud = —2,P the Hamiltonian reads 



H2 = J'[(5l, - ^2.)' + {Sly - S2yf - S2 - S2]. (12) 

Applying the unitary (rotation) transformation 
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FIG. 2: (Color online) Ground state probability density 



p(nd) = IM^dW for + 1 = 10"^ (a) and 



+ 1 = -10" 



(b). Solid black lines are related to the number of particles in 
each gas component 2N = 50, dashed red lines to 2N — 100 
and dotted-dashed blue lines to 2N = 200. 



energy manifold (n^ = 0). That is, all values of Ud are 
equally probable. 

It is interesting to note that we can construct the exact 
quantum ground state using the superposition of sym- 
metry broken solutions obtained in the classical limit. 
At the transition point the classical analysis tells us 
that in the ground state the potential wells are associ- 
ated with classical angular momenta where orientation 
of one is given by {9i, (pi) = {6, cj)) and the other one by 
(6*2, <t>2) = (■"■ — ^, 0) and spherical angles 9 and 4> can be 
arbitrary. The best quantum approximation of a classical 
angular momentum is a coherent state psl [26| 



which commutes with 5j2 and thus leaves the manifold 
5i2 + 522 = invariant, and defining the total spin oper- 
ator, S = Si -|- S2, we can rewrite the Hamiltonian ()12|) 
in the following form 



i?2 = J'(S2-S2-S2) 



(14) 



States |j, 0) with an integer j, where j{j + 1) is an eigen- 
value of the operator and is an eigenvalue of the 
Sz = 512 + 5*22 operator, are therefore eigenstates of our 
problem. The energy spectrum reads 



E,=j\j{j + l)-2N{N + l)), 



(15) 



with < i < 2N and the ground state solution corre- 
sponds to i = 0. In the basis ^ it takes the form of 

1 ^ 

|^o> = -/^^= \N+n,N-n)\N~n,N+n). (16) 

n—~N 

Huge particle number fiuctuations become apparent in 
Eq. (|16p where the ground state turns out to be a uniform 
superposition of all Fock states belonging to the lowest 
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(17) 



If we postulate that the quantum ground state of our sys- 
tem can be approximated by a single tensor product state 
|0, 4))\'K — 9, (j)) the rotational symmetry of the Hamilto- 
nian (jlOp will be broken. A rotationally invariant state 
can be restored if we prepare a uniform superposition of 
the tensor product states, i.e. by integrating over all solid 
angles 



d(f> 
2^ 



sin9d9\9,(f>)\Tr-9,(t)), 



(18) 

where (V'olV'o) = 1- Substituting Eq. dTT]) into Eq. (fT8|) . 
integrating over </> and employing the identity 
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FIG. 3: (Color online) Ground state probability densities 
for a number of particles in each component 2N = 1000, 
J = 5 X IQ-^ and ^ = -0.9999 (a), ^ = -1 (b) 
and -^j^ = —1.0001 (c). Solid black lines correspond to 
Pi'n-d) = X^n^ \'^o{ns,nd)\^ where ^po{ns,nd) is obtained by 
diagonalization of (O while dashed red lines are related to 
p(nd) = \ipoind)'^ with ipo{nd) obtained by diagonalization of 
the effective Hamiltonian (|10p . In panels (a) and (c) the black 
and red lines are hardly distinguishable. 

that follows from the completeness relation of the coher- 
ent states we restore the quantum ground state (fT6|) . 

We see that at the transition point the classical anal- 
ysis allows us to construct the exact ground state of the 
system. However, in the close vicinity of the transition 
point this analysis is not able to provide a good esti- 
mate for the ground state of the system. This is because 
the classical approach predicts discontinuous transition 
between the phase separation and miscible regimes but 
the transition is actually continuous. Starting with the 
classical ground states and using the coherent states we 
can construct quantum ground states that depend on N 
and on the sign of the parameter Ud/{2J'^) + 1 but not 
on its absolute value. However, the exact diagonaliza- 
tion indicates that there is a range of Ud/ (2 J^) -I- 1 where 
the ground state changes continuously from the misci- 
ble to phase separation character. As one can expect 
this range shrinks with N because differences between 
the classical and quantum angular momentum diminish 
in the large N limit. This is illustrated in Fig. [2] where 
wc plot p{nd) = \Mnd)\^ for Ud/{2J^) + 1 = lO'^ and 
— 10~^, i.e. close to the transition point, for different 
values of N obtained in numerical diagonalization of the 
effective Hamiltonian ((TO)) . 

D. Validity of the perturbation approach 

Our predictions arc based on the effective Hamiltonian 
(llOp which is second order in J, and they are correct pro- 



vided higher order terms can be neglected. The fourth 
order terms are smaller than j4(27V)'* or j'^\Ud\{2NY. 
At the transition point \Ud\ = so if J'^{2NY is 

much smaller than the energy gap between the ground 
and first excited states of the Hamiltonian ([T4)). i.e. 
El — Eq = 2J-^, then the higher order terms can be 
neglected and the system is properly described by the 
second order Hamiltonian ([Ti)) . Hence, 

J^{2Nf < 1, (20) 

is a sufficient condition for the validity of the ground state 
that describes the huge particle number fluctuations 
in the system. 

Wc have tested our predictions comparing them with 
exact numerical calculations. Figure [3] shows probability 
densities p{nd) corresponding to ground states of the sys- 
tem obtained by diagonalization of the full Hamiltonian 
^ and the effective Hamiltonian ^ for 2A^ = 1000 and 
J = 5 X 10^^. Different panels are related to different 
values of Ud in the vicinity of the transition point. That 
is, -f 1 = 10""* corresponds to the miscible regime, 
+ 1 = to the transition point and ^ + 1 = -10"'' 
to the phase separation side of the transition point. We 
can see that on both sides of the transition point the 
density p{nd) is peaked around the classical solutions, 
while at the transition point it is uniformly distributed. 
Figure [3] indicates perfect agreement between the per- 
turbation calculations and the exact results even though 
J^{2N)^ = 0.25 and thus the condition (gHl) is barely 
fulfilled. 



IV. CONCLUSION 

In summary, we have analyzed a strongly interacting 
two-component Bose gas in a double well potential for 
parameters close to the transition point where the phase 
separation occurs. The second order effective Hamilto- 
nian allows us to describe the system in the vicinity of 
the transition point when higher order corrections are 
negligible. We have shown that, at the transition point, 
the ground state of the system becomes a uniform su- 
perposition of Fock states. That is, the system reveals 
huge quantum fiuctuations of populations of the poten- 
tial wells. 
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